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ABSTRACT 



This thesis explores the use of a pointed seeker lens 
designed using a spherically symmetric gradient refractive 
index (GRIN) . The design helps to solve the current design 
conflict between optical quality and aerodynamic drag in- 
herent in hemispherical seeker lenses. Equations for lens 
design and the evaluation of off-axis lens performance have 
been developed for both a homogeneous version and a GRIN 
version of the pointed seeker lens. The homogeneous lens 
is used as a comparison and a check for the GRIN lens. A 
FORTRAN program (GISL) has been written and employed to 
evaluate and compare both the homogeneous lens and many dif- 
ferent configurations of possible GRIN lens designs. Results 
indicate that the GRIN lens has highly superior off-axis imag- 
ing performance as compared to the homogeneous lens. The best 
results were obtained for the GRIN lens with a fifty percent, 
positive, spherically symmetric gradient index with center of 
symmetry interior to the lens. Only very slightly inferior 
performance was observed with a five percent version of the 
same lens; such a lens possibly can be manufactured today. 

GRIN lens performance also indicates that for objects off- 
axis by more than 17.2 degrees a large scale, multiple ele- 
ment sensor array may be required; with such a sensor array, 
objects off-axis by more than 37.2 degrees may require mirror 
elements to compensate for image movement. 
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I. INTRODUCTION 



Historically, the design of tactical missiles employing 
passive or semiactive infrared (IR) seekers has involved a 
difficult compromise between aerodynamic requirements and opti- 
cal or seeker requirements. Whereas aerodynamically the mis- 
sile nose region should be sharp in order to reduce drag, 
optically it should be hemispherically shaped for image 
quality and as large as possible to increase aperture and 
therefore acquisition or tracking range. Some IR homing mis- 
siles designed for very short range anti-armor missions have 
totally ignored nose drag in order to optimize seeker per- 
formance while other designs for longer range missiles re- 
quiring high cruise velocities and greater aerodynamic 
efficiency have used the blunted ogive as a compromise. There 
have not been any IR designs which have ignored optics in favor 
of aerodynamics; nor has there been employed a pointed seeker 
lens with the desired optical qualities. 

In order to increase the performance of optically guided 
missiles beyond the current state of development, the conflict 
between aerodynamic requirements and optical restrictions must 
be resolved. Significant improvement in missile performance 
by increasing thrust is not likely due to the highly advanced 
state of propulsion today. One way to resolve the aerodynamics- 
optics problem is to design a pointed lens which has, if not 
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imaging quality, enough optical performance to allow the 
reduction of tracking data. 

The lens must be in the general shape of a cone or ogive 
with half angle compatible with the design free stream Mach 
number in order to maintain an attached shock wave. Such a 
spike— shaped lens might also be used in the diffuser portion 
of an optically guided ramjet with nose inlet to conserve 
stagnation pressure. A diffuser- lens may benefit by a semi- 
isen tropic spike shaped lens. 

Poor optical performance due to the pointed shape of 
sharp lenses has precluded their use. Gradient Refractive 
Index (GRIN) materials, however, permit the lens designer 
the freedom to spatially vary the lens index of refraction 
to compensate for a traditionally poor optical shape. Al- 
though the use of GRIN has seen widespread use in fiber optics 
technology, it has been used infrequently in lens applications 
until recently. Lens designers are discovering that multiple- 
element photographic objectives may be redesigned using a 
two-element gradient lens [1]. It should be noted, however, 
that such GRIN lens have not yet been successfully fabricated 
even though large index changes in glass have been accomplished 
by the diffusion of doped electropolarizable ions [2]. Ex- 
tensive research is being conducted in the creation of ever 
larger and more precise gradients. A total change in refrac- 
tive index of approximately five percent is currently 
attainable . 



37 



At the Naval Postgraduate School, Frazier [3], Terrell 
[4], and Amichai [5] have studied GRIN as applied to the 
sharp lens problem. The brief introduction by Frazier was 
followed by Terrell who designed a sharp lens having a coni- 
cal outside surface and a variable inner surface using a 
homogeneous index (HIN) and then also briefly touched on the 
GRIN application. Amichai extended Terrell's lens to the 
GRIN case and included preliminary results from this thesis 
in a computer routine intended as a framework for lens 
optimization by following researchers. 

This thesis investigates the design and off-axis per- 
formance (skew rays) of a GRIN seeker lens having a variable 
outer surface and a fixed right circular cone as an inside 
surface. The variable outside surface is determined by the 
character of the spherical gradient employed and varies from 
a pseudo ogive to an approximate isentropic spike as shown in 
Figure 1. First, the theory of the HIN lens is developed, 
followed by the design of a homogeneous lens intended both 
as a comparison and a check for the GRIN lens theory and the 
design which follows. Lens performance parameters are dis- 
cussed and results presented for both the HIN and GRIN lenses. 
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OUTSIDE CONE 




(a) Terrell's Lens With Conical Outer Surface 




(b) Lens Shapes In This Thesis With Conical 
Inner Surface. 



Figure 1, Lens Shapes 



39 



* . 
I.j 



II. THE HOMOGENEOUS LENS 



A. THEORY 

Snell's law is the cornerstone of contemporary lens de- 
sign in homogeneous optical materials. In the HIN case 
Snell's law is used in the familiar form 

sin sin 0^ (1) 

where is the index of refraction of the material corres- 
ponding to surface intercept angle 0^, with respect to the 
surface normal, at the interface between surfaces. Lens 
geometry and the relative values of N^ determine the resultant 
optical behavior. 

B. ASSUMPTIONS AND SIGN CONVENTION 

In order to simplify the design and analysis of the 
seeker lens problem, certain assumptions have been intro- 
duced. Although energy loss upon transmission through the 
lens at each surface is calculated, it is assumed that the 
light is monochromatic radiation, time dependent 
electric and magnetic fields. Light impinging upon the lens 
is assumed to have a planar wave front as if propagating from 
an object at infinity. The presence of a shock wave attached 
to the lens is ignored as are any other regions of expansion 
or compression in the flow field about the lens [6] . Further- 
more, the index of refraction of the free stream is assumed 
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to be equal to that of the interior space behind the lens 
even though computer routines were written with the flexi- 
bility to process unequal values . 

The sign convention used is a right handed system with 
spatial coordinates positive to the right, up, and out of 
the page as seen by the reader. All angles are assumed 
positive counterclockwise from point of reference unless 
otherwise noted. All linear dimensions are implicitly non- 
dimensionalized with respect to the maximum radius of the 
right circular cone forming the inside surface of the lens. 

C. HIN LENS DESIGN 

The lens design procedure consists of calculating both 
the loci of points forming the outside surface and the slope 
of the surface at each point in the meridian plane. It is 
convenient to approach the problem by placing a point source 
of light at the design focal point F and calculating succes- 
sive refracted ray paths U' using Snell's law at points 
on the inside surface in Figure 2. Points E^ are formed by 
the intercept of the refracted ray and the slope of the out- 
side surface as extended from the previous point The 

tangent to the surface which will refract the ray in the de- 
sired direction, U, may again be found from Snell's law. 

The accuracy of the calculated surface increases as the 
spacing between points decreases or as angle AU" becomes very 
small. The final ray is parallel to the lens axis and point 
corresponds to point B, which is the origin. The area 
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Surface Solution 
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formed by the triangle ABE^ is opaque to prevent rays from 
scattering through the opposite side of the lens. 

A particular lens is begun by specifying the focal 
length BF , the cone half-angle a, the indices of refraction 
f and the cone radius at the edge R, the ray 

direction angle U, and the total number of rays to be traced. 
The thickness, T, of the lens at the edge (E^D^) must also 
be specified in order to define which of the family of possi- 
ble outer surfaces will be calculated. In general, angle U 
will be taken to be zero in order to investigate objects on- 
axis at infinity. R will always be set equal to one and to 
nondimensionalize all linear dimensions every length is im- 
plicitly expressed as a ratio with respect to R. The angle 
3 at the focal point is measured from the lens axis to and 
is expressed as 



tan I 6 



R 

(BF - R'cotan a) 



• ( 2 ) 



Thus, if I+l is the total number of rays to be traced, then 



AU" = S/I 



( 3 ) 



and 



U" = 6 - J‘AU" 



( 4 ) 



where J is the ray number. To begin the lens design the 
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coordiantes X2 ( J) lY 2 point Dj must be found. Here 
the subscript 2 refers to the inside surface where x^(J) , 
Vj^(J) are the coordinates of Ej on the outside surface. 

rollowing Kingslake [ 7 ] , Q in Figure 3 may be expressed 
as 



Q = BF sin U" = X2(J) sin U” + 72 (J) cos U" ( 5 ) 



But , 



72 (J) ^ 



( 6 ) 



so that 



BF sin U" = X2(J) [sin U" + tan a cos U"] ( 7 ) 



Now, X2(J) and 72 (J) are 






BF sin 0 " 

[sin U" + tan a cos U”] 



( 8 ) 



72 (J) 



BF in U" tan a 
[sin U" + tan a cos U"] 



( 9 ) 



Looking now to find the incident angle which the ray forms 
with the inside surface at Dj it can be seen that 



tan (a + I2) 



BF - X2(J) 



( 10 ) 
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Figure 3. Ray Trace Geometry 



Substituting into Equation (1), Snell's law, the refraction 



angle can now be written as 



I 



2 




( 12 ) 



The path of the ray inside the lens may now be expressed 
in the form of the angle U' where 



To this stage all rays are treated the same. Now, however, 
a differentiation must be made between the outermost ray 
which defines the edge of the lens and all other subsequent 
rays . 

For the outermost ray the thickness, T, of the lens at 
the edge can be specified which yields x^(l) and y^(l) 
immediately : 



U 



(2 - a) - I 2 



(13) 




X 2 (l) - T cos U’ 



(14) 



y^(l) = y 2 (l) + T cos U' 



(15) 



At point we may write Equation (1) again to obtain 



sin ^ [ ( jy^) sin 



(16) 
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and are the incident and refracted angles 
at the outside surface and are not yet known. Angles 
and may be found by first noting that 



U + = U' + 



(17) 



then substituting Equation (16) into (17) and rearranging 
terms to obtain 



-1 



N. 



= U' - U + sin [(^)sin 



(18) 



It can be shown after some algebraic manipulation that 
Equation (18) may be solved for the angle in the form 



= 3in“^{ 



sin (U' - U) 



[cos (U'-U) - N^/N 2 l ^ + sin^ (U' - U) 



1/2 

■} (19) 



Angle may be found, if desired, by substitution of Equa- 
tion (19) into Equation (16). Of more importance, however, 
is the determination of the slope of the outer surface and 
the surface normal. Now that angle is known these slopes 
may be written as 



= 

dx I normal 



DYDXN(J) 



tan(U + Ij_) 



( 20 ) 



DYDXT ( J) 



dx 



surface 

tangent 



cotan (U + 



( 21 ) 



Equations (20) and (21) may be evaluated by using Equation 
(19) for 

Now that all the parameters are known for the first ray, 
the remainder of the points E^ and the respective slopes 
may be found. Each successive ray is traced as before by 
Equations (4) through (13) . Equations (14) and (15) , how- 
ever, may not now be used since the lens thickness along the 
ray can not be specified. Instead, the intersection of the 
ray and the surface slope from the previous ray is used to 
define the new point The intersection is found by first 

writing the equations of lines representing the ray and the 
surface tangent. For the ray: 



y = -X tan U' + 72^^) tan(U') (22) 

For the surface tangent; 

y = X cot(I^+U) + y^(J-l) - Xj^(J-l) cot(I^ + U) 

(23) 

Equations (22) and (23) are solved simultaneously to yield 
the coordinates of E^ which are 

yj^(J) = (24) 
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r 





where ; 



A = cot(I^+U) cot(U') tY 2 (J) +x^ (J) tan(U' ) ] (25) 

B = y^(J-l) - x^(J-l) cot(I^ + U) (26) 

C = 1 + cot(I^ + U) cot(U') (27) 

and 



Xj^(J) = cot (U') [-Yj_ (J) +Y 2 (J) +X 2 (J) tan(U' ) ] (28) 



Now Equations (19) , (20) , and (21) are used to calculate 
and the slopes at Ej. Therefore, the remainder of the lens 
surface may be generated. The opaque region at the surface 
is formed by extending the slope of the surface at E to 
intersect the lens axis. Here, K refers to the last ray. 
The nose half-angle, y, thus formed is given by 



tan y 



= 

dx I surface , 



K 



(29) 



D. SKEW RAYS 

A skew ray is one that begins from an off-axis object 
point and enters the lens either in front of or behind the 
meridian plane (z = 0) . For every skew ray, there is a 
corresponding mirror image skew ray on the opposite side of 
the meridian plane so that two skew rays are traced at the 
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expense of only one calculation. These two skew rays inter- 
sect at the same diapoint. 

Large numbers of skew rays are traced through a lens 
in order to study lens performance at different obliquities. 
The procedure is to superimpose a grid over the lens aperture 
and to trace rays through the intersections of the grid, 
through the lens and onto the image plane. An image plane 
spot diagram and an energy density plot are then constructed 
for study. 

To accomplish this, the aperture grid has been attached 
to the nose (opaque region) of the seeker lens at station A; 
see Figure 3. The plane of the aperture grid is tilted 
relative to the lens by a variable angle. A transformation 
between the grid coordinates and the lens coordinates has 
been derived to connect skew rays from grid to intercept with 
the outside surface of the lens. Referring to Figure 4, it 
can easily be seen that 



X 



I 



y 



1 



z 



» 



[x + AB]cos a - y 
P 

[x + ABjsin a + y 

iT 



Z 



sin a 

P 



cos a 

P 



(30) 

(31) 



(32) 



where a is the tilt of the grid plane, AB is the length of 
P 

the opaque nose portion on the x-axis. Both z and z' are 
positive out of the page. 
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Figure 4. Coordinate Transformation (z and z' are out of the page) 



The three dimensional lens outer surface is generated by 
rotating the array of outer surface coordinates in the merid- 
ian plane through 2u about the lens axis. Each pair of 
coordinates x^(J),y^(j) thus describe a circle in the y-z 
plane of the lens. This circle transforms, however, to an 
ellipse in the grid plane given by the equation 



y ' 



[x + AB] 



sin a 

P 





cos 



a 



P 



(33) 



where R is the radius of the circle in the x-y plane. R 
may be expressed by the familiar equation of a circle 



R 



2 



2^2 
y + z 



(34) 



in the y-z plane of the lens, or by solving Equation (33) 
for R in the grid plane 

= z'^ + ( — - (x + AB) tan ot ) ^ (35) 

cos a p 

P 

With the aid of the foregoing groundwork, the x,y,z coor- 
dinates of the ray intercept with the outer surface may be 
found. To see how this is accomplished, first refer to Figure 
5. The skew ray will pass outside the circle formed by rotat- 
ing some point x^(J) ,yQ^(J) and inside the next circle formed 
by rotati.ng x^ ( J-1) ,y^_ ( J-D . In so doing the skew ray will 
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DYnXN(J-l) 



X 
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Figure 5. Geometry of Skew Ray In tercept_Wi th Outside Surrace. ^ 

Plane is the Meridian Plane. is the Surface Normal in the 

Meridian Plane; N' the Surface Normal at the Point of Intercept 

PI. 



intercept an imaginary cone formed by the two circles.^ 

Now, the y',z' of the skew ray in the grid plane and a 
trial value ( J) from the array of surface points are sub- 

stituted into Equation (35) and a value for R is calculated. 
If R is greater than yj_(J) but less than y^(J-l) than the 
appropriate circles have been found. If not, another trial 
value of Xj^(J) is picked. In practice the x^(J) corresponding 
to 



y3^(J) = 2' (36) 

is chosen as the first trial value and subsequent trial 
values are x^(J-l), x^(J-2) and so forth. 

The equation of the cone passing through the two circles 
may be derived by again referring to Figure 5 and noting 
that 



QP = Qp' = v'y+z = R (37a) 



QP = y = m[x - x^(J)] + y^^(J) 



(37b) 



where the slope m is given by 



^Recall that a straight line between points Ej was used 
to form the outer surface in the meridian plane. These 
approximations are valid only if the spacing between points 
or circles is very small. 
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m 



(38) 



y3^(J-l) - y^(J) 

x^(j-l) - Xj_(j) 



and therefore 

2 y-, ( J-1) - y, (J) 2 

3y equating Equation (39) to Equation (35) and solving for 
X, the expression for the x-coordinate , Xq , of the outside 
surface ray intercept may be found. After some algebra, this 
reduces to the complicated relationship 



PARI ^ ypARl _ PAR3 
^0 PAR2 " PAR2 PAR2 



where : 

PARI = ab + ce - b^AB - c^d (41) 

PAR2 = b^ - c^ (42) 

PAR3 = a^ + b^ (AB) ^ - 2abAB + z ' - (e-cd) ^ (43) 

and where 

a = y'/Cos ot 

P 

b = tan (4 5) 
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c = 


1 1 

1 1 

X k; 


d = 


x^(J) (47) 


e = 


Yl(J) (48) 



It has been found that the plus sign in Equation (40) gives 
the correct values until PAR1/PAR2 becomes greater than 
Xj^(l) at which time the negative sign must be used. It now 
follows from Equations (31), (32) and (40) that the and 

Zq coordinates of the intercept are 



o 

II 


■ (^0 “p 

P 


and 





N 

O 

II 


z’ (50) 



The next step in tracing the skew ray is to ascertain the 
direction cosines of the ray inside the lens after refrac- 
tion at the outside surface. This is accomplished in a 
series of steps beginning with the determination of the 
angle ^ which the skew ray makes with the normal to the 
surface at the point of intercept. Angle <p rnay be found by 
taking the scalar product of the direction cosines of the 
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skew ray and the surface normal. Thus to find (p , the surface 
normal must first be obtained. 

As might be sux'mised by the reader, the direction cosines 
of the surface normal can be found by taking the gradient 
of Equation (39) with Equation (37) substituted for R. 

Several orders of magnitude may be gained in accuracy, how- 
ever, if the normal is found by interpolation in the 
meridian plane and then rotated to the point of intercept by 
angle u as in Figure 5. This interpolation procedure involves 
multiplying the ratio S/ST with the difference between 
DYDXN(J) and DYDXN(J-l) and adding the product to DYDXN(J) . 

DYDXNpjj^ = ^(DYDXN(J-l) - DYDXN( J) ) + DYDXN( J) (51) 

where subscript PIM refers to point of intersection, meridian. 
Values of the slopes of the normals in the meridian plane, 
DYDXN(J) , are given by Equation (20) and S and ST are the 
linear separation of points P", PIM and P", P respectively. 
S/ST is given by 



Therefore, the normal vector, N^, in the meridian plane is 



[(R-y^(J))^ + (XpiM “ ^ ^ 




S 

ST 




(52) 



N. 



m 




(53) 
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and since 



tan \i 



z ' 



^0 




(54) 



the normal vector at the point of intersection on the 
outside surface becomes 



N' = ( 



n vnxN cos (tan j + sin (tan k (55) 

PIM ^0 



^o' 



and 



o 1/2 

|N’l = [DYDXNpjj^'^ + 1] (56) 

Now the scalar product of the skew ray vector, R, and surface 
normal, N', may be performed to find (}) . Since the acute 
angle between these vectors is required, the dot product 
must be written 



R-N’ = |N'l cos(tt-'^)) 

Substituting Equation (55) for N' and noting that 



R = Ki + Lj + Mk = cos ttpi - sin a^j 



the left side of Equation (56) may be expanded to 



• M ' = 



R-N 



COS a .1 Zq 

DYd 3^- sinapCos(tan -) 



(57) 



(58) 



(59) 



58 



Solving Equation (56) for the angle <J) and introducing 
Equation (5 8) / the expression for angle <() becomes 



cos a 



(}) = TT-COS 



-1 



DYDXN 



^ _ 



PIM 



sin a cos (tan 



-1 ^ 0 , 



[DYDXNpjj^‘2 + cos^(tan'^ + sin "( 



I 

in- 2 ( tan- 1 !i) 1 1/2 I 



/O 



(60) 



Now that the acute angle, cj), between the ray and the outside 
surface normal is known, (}> ' , the angle between the refracted 
ray ihside the lens and the surface normal may be found using 
Snell ' s law. 



<!>' = sin"^(^ sin 4>) (61) 

2 

The next step is to ascertain the direction cosines K', L', 

M' of the refracted ray inside the lens. Following Kingslake 
[7] , Figure 6 shows the optical vector relationship between 
R, the skew ray, R' , the refracted skew ray, N', the surface 
normal, and the indices of refraction of the two media. 
Algebraically this relationship is written as 

^2^' ^ (N 2 cos cos (j)) N' (62) 

By resolving Equation (61) into component form, the direction 
cosines of R' may be found. After rearranging, K', L', and 
M ' be come 
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Figure 6. Kings lake's Skew Ray Diffraction 
Diagram 
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N. 



I 1 

= N~ ^ cos (}))k 

2 



( 63 ) 



^1 ^1 
= N~ ^ '^' ■ FP ^ 



(64) 



M' 



Ni N 

= n" M + (cos <{)'-_ cos 4)) 



m 



(65) 



where K, L, and M are given by Equation (57) and k, m are 



k 



Z 



In' 1 



cos (tan ^ — ) 
^0 




( 66 ) 



(67) 






sin ( tan 



-1 





( 68 ) 



Again, |N | in Equations (66) through (68) is expressed by 
Equation ( 56) . 

At this stage in the process of tracing the skew ray, 
the coordinates of the external ray intercept Xq , y^ , 
are known as are the direction cosines K', L', M' of the 
skew ray inside the lens. The next step is, of course, to 
find the ray intercept with the inside surface (cone) and 
the direction of the ray subsequent to refraction. The 
intermediate steps are similar to, if not identical with. 
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the foregoing. One major difference, however, is that the 
inside conical surface may be expressed analytically and the 
intercept coordinates y^, may be found exactly without 
approximation. Proceeding, it can be seen that if x. , v., 

2 ^ were already known. Equations (63) , (64) , and (65) could 
be rewritten as 



K’ 




(69) 



L’ 




(70) 



iM' = 




(71) 



where D is the distance between surface intercept points. 
Upon rearranging: 



Xj_ = DK' + Xq (72) 

y^ = DL' + yQ (73) 

z.=DM'+z- (74) 

1 0 

Furthermore, there exists a relationship between x^, y^^, and 
2 ^ which is given by the expression for the inside conical 
surface ; 



y^ + z? - X? tan^ a = 0 (75) 

^ 1 . 1 1 
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where ex is the half-angle of the cone. Substituting Equa- 
tions { 12 ) , (73) , and (74) into Equation (75) and solving 

for D, it is seen that 



D 





- P P 
2 3 



(76) 



where 






2 

- K'Xq tan a 



(77) 



?2 = L’2 + m’2 - k' 2 tan^ a 



(78) 



and 



= Yr 



2 

"o - 



2 ^ 2 
Xq tan a 



(79) 



Correct values for D are obtained by using the minus sign in 
Equation (76) . Now that D is )cnown, values for K', L', M', 

Xq , Yq , and Zq are substituted into Equations (72) , (73) , and 
(74) to yield the coordinates of the inside surface intercept 
point. 

In the case of the inside surface, the gradient may be 
used to obtain the surface normal at the point of intercept 
since this surface has been expressed analytically. If 
Equation (75) is denoted by f(x,y,z) then 



N. 



1 



7f 



7f 



X . 
1 



■Yi. 



z . 
1 



(80) 
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where is the unit vector in the direction of the surface 
normal at the point of intercept. Here, 



9f;^8f^^9fr 
9x 9y 9z 



or alternately. 



7f = fi + fj+fk 
X y-* z 



(81) 



(82) 



Thus, the direction cosines of the normal k', I', m' may 
be expressed as 



^ 

Vf I 



m 



Vf I 
f 

z 

Vf 1 



where : 



= - 2Xj^ tan a 



f = 2y ■ 

y ^ 



fz = 2z. 



(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 
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and 



1 — , 2 4 2 2 ^ 

|Vf| = 2 [x^ tan cx + + zf] (89) 

^ /V 

Now that the two unit vectors R' and N' are known the inci- 
dent acute angle of the skew ray with the inside surface, 

may be ascertained by again taking the scalar product as 
in Equation (57) . Thus, 

= TT-cos”^ (K'k' + L'i' + M'm'] (90) 



and using Equation (1) again 



-1 ^2 

= sin tf-sin ^ 



(91) 



where is the acute angle between the skew ray and the 
surface normal after refraction at the inside surface. 
Direction cosines K", L", M" of the skew ray after refraction 



are 


found 


analogously to K', L' 


, and M' 


' in 


Equations 


(63) , 


(64) 


and 


(65). Here, however, 


N 3 , 


K', 


L ' , M' , 4> . 

1 




and 


k' , i 


' , m' are substituted 


for , 


N 2 , 


K, L, M, 




and 


k, 1, 


m respectively. Thus 











K" 



L” 



^2 ^2 

_ K- + (cos ^ cos 4>i)k' 

^2 ^2 

— L' + (cos 4>: - ^ cos 



(92) 



( 93 ) 
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and 



M" 



2 

^ M ' + { cos <J) 



cos (j) . ) m’ 
1 



(94) 



Finally, the skew ray intercept with the image plane may 
be found. The image plane is treated as another surface 
along the path of the skew ray, and the intercept coordinates 
^im' ^im' ^im easy to find. Equations (72) , (73) and 

(74) may be used again in the form 



X 



im 



D'K" + 



(95) 



im 



= D'L" + y. 



(96) 



z . = D'M" + z . (9 7) 

im 1 



where D', in this instance, is the linear separation between 

coordinates x^, y^, at the inside surface and , y^^, 

z. of the image olane. Furthermore, 
im 



X 



im 



BF 



(98) 



is the equation of the image plane. It follows that 



D' 



(BF - x^) 



(99) 
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after substitution of Equation (98) into Equation (95) . 
Therefore , 

BF - X. 

and 

BF - X. 

"im ' > — noi) 

Hence, the skew ray has been traced onto the image plane 
at the focal point of the lens. The y. coordinate of the 
corresponding mirror image skew ray is the same as Equation 
(100) . The coordinate, however, is the negative of 

Equation (101) since the mirror image ray is behind (when 
viewed along the z-axis) the meridian plane. After tracing 
a complete set of skew rays through the lens a spot diagram 
may be plotted. Clearly, the number of rays to be traced 
depends entirely upon the incremental size of the aperture 
grid chosen. According to Kingslake [7] , at least 100 rays 
must be traced to give a fair approximation of the actual 
image. In addition to the spot diagram, an energy density 
plot may now be constructed by counting the number of rays 
within progressively larger radii from the image centroid and 
then plotting the number of rays as a function of radius. 
Here, each ray is assumed to contain a unit, nondimensional , 
an amount of radiant energy for convenience. 
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E. RADIANT ENERGY LOSS 

In reality , each ray loses intensity upon transmission 
at each interface. Of the total amount of energy contained 
in each ray, a fraction I will be transmitted, a fraction I 

R 

will be reflected, and a fraction I^ will be absorbed by the 
medium into which the ray is propagating. Since it has been 
assumed that absorption is negligible, it must be true that 



Furthermore, the relative amounts of transmitted and reflected 
electromagnetic energy may be calculated by the well-known 
Fresnel Equations which state the dependency of I^ and Ij^ 
upon the angle of incidence and the indices of refraction 
at the interface. That I^p and 1^^ are further dependent upon 
the orientation of the electric vector with respect to the 
geometry of ray incidence is fundamental to the boundary 
conditions which govern the form of the Fresnel relations as 
derived in Hecht-Zajac [8] . Since this thesis examines lens 
response to monochromatic radiation, the indices of refraction 
are not considered as a function of wavelength; further, the 
electric vector orientation is assumed to be rapidly and 
randomly changing with time. By time averaging field com- 
ponents, it may be seen that the reflectance is 



I, 



T 



+ I 



R 



1 



( 102 ) 



I 




( 103 ) 
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where 



r 




(104) 



and 




(105) 



Here, is the ratio of the index of refraction of the 

transmission side of the interface to the index of refrac- 
tion of the incident side. From Equation (102) it now follows 
that the transmittance through the interface is 



The total transmittance through the lens is simply the product 
of at the outside surface with that of the inside surface 
where is substituted for g in Equations (104) and (105) . 

Total internal reflection of the ray may occur at the 
inside surface if the incident angle becomes too large. 
Following Reference (7), this occurs when angle is equal 
to or greater than it/ 2. Thus, Snell's law becomes 



I 



t 



1 - I 



R 



(106) 



N 



3 



sin g . 



N 



(107) 



2 



and any ray with equal to or greater than this will be 
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totally internally reflected. In this thesis such rays are 
labeled "failed rays" since they fail to intersect the image 
plane. If N 2 is 1.5 and is 1.0, the incident angle for 
total internal reflection is 41.81°, or greater. 

F. OPTICAL PATH LENGTH (OPL) 

The optical path length of a skew ray is an analytical 
tool with which the researcher may ascertain the phase of a 
ray at the end of the path. By so doing, the image diffrac- 
tion pattern may be constructed which shows the addition or 
subtraction of amplitude depending upon relative phase. Since 
each ray must begin with the same phase, monochromatic radia- 
tion is used for diffraction experiments. OPL is included 
here only as a matter of interest. The calculation of optical 
path length is simply the sum of the geometrical path segments 
of a ray multiplied by the corresponding index of refraction 
of the medium, for that segment. 
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III. THE GRIN LENS 



A. THEORY 

This thesis assumes a spherically symmetric, inhomogene- 
ous, isentropic medium in which the refractive index varies 
from point to point but is independent of direction at each 
point. The refractive index is a function of the coordinates 
of the points of the region being considered. The problem 
of describing the resulting curved paths of rays in such 
GRIN materials has been solved long ago in the form of a 
single second order vectorial differential equation. Marchand 
[9] has shown that the solution to the differential equation 
in the case of spherical gradients can be written in polar 
coordinates in the plane of the ray as 



9 



+ e 



dr 



r 2 2 

r (n r 



.e2]l/2 



( 108 ) 



Here r^ and 9 q are values of r and 9 at a convenient reference 
point on the ray and ]e[ is a scalar constant along the ray 
given by 



e = £ n^ r^ sin (109) 

Referring to Figure 7, r is measured from the center of 
symmetry of the index function; angles 9 and 9 q are measured 
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Figure 7. Grin Lens Geometry 
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are 



counterclockwise from the lens axis; angles ip and 
measured counterclockwise from the radius vector to the 
instantaneous ray direction vector. £ is the sign function 
expressed as 

£ = sgn = ±1 ( 110 ) 



where 



I’J'ol < ^/2 (111) 

In Equation (109) , n^ refers to the index function n evaluated 
at r^ . The form of n has been chosen as 

n = n(Og,r) (112) 

where denotes the position of the center of symmetry on 
the lens axis (shown in Figure 7) and r is the radial coor- 
dinate from 0^ to the point in question. This form allows 
the study of the effects on lens performance as 0^ is changed. 

More specifically, in order to allow the analytical 
integration of Equation (108) , the expression used for 
n(0^,r) is a generalized version of that employed by Luneburg 
as described by Marchand [9]. Here 

2 

n = [a + b(^) (113) 
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II 



whereas Luneburg used specific values for a and b. Generaliz- 
ing the index function enables the strength, or percent change 
of the gradient to be varied as well as the algebraic sign 
of the gradient. If the parameter b is negative, a decreasing 
parabolic gradient results; conversely a positive b yields 
an increasing parabolic gradient. By visualizing a plane 
wave front passing through a spherical gradient, it may be 
seen that the negative gradient results in light rays bending 
forward at the center, as shown in Figure 7, whereas a posi- 
tive gradient has the opposite effect. This principle was 
used by Wood, as related by Marchand [9], in constructing 
simple lenses having plane faces and a radial index. The 
wood lens acted as a converging or diverging lens depending 
upon the sign of the gradient used. 

Additionally, it must be noted that if the parameter b 
is equal to zero. Equation (113) reduces to 

n = /I (114) 

which is a HIN lens having constant index of refraction. 

This fact has facilitated the correlation of GRIN and HIN 
computer trace algorithms. 

Equation (108) may be integrated by a change of variable 
using the relation 
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*■ 



which leads to the solution 



9 = 



9 - — 

2 



sin 



-1 



2e^/r^ - a 
^/a^ + 4be^/R^ 



o 2 , 2 

-1 r 2e /r. - a 



- Sin 



va+^be/^-* 



Z 

(116) 



Equation (116) gives 9 as a function of r. This equation 
can be easily solved for r as a function of 9 in the form 



r = 



a+- 



/2 



I 2 ^ 4be . 
a H ^r-sin 



r: 



-2e (9-9g ) +sin 



-1 



o 2 , 2 

2e /r. - a 



V^a^+4be^/R2 



1/2 



(117) 

The instantaneous direction of the ray at any point r,9 
may be ascertained by using the invariance of e. Hence, 



e = 



;nrsin’(; = en^rgSin'^ 



(118) 



and 



'P 



. -l,"o "^0 *0, 

sin [ zn:: J 



nr 



(119) 



Furtherirore , the orientation of the plane of the ray may be 
easily deduced since every ray in a spherical medium is a 
olane curve lying in a plane through the center of symmetry. 
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Using this fact, Marchand [9] has shown that a suitable 
conversion from coordinates r and 6 in the plane of the ray 
to global Cartesian coordinates may be written in the form 



X 



y 



z 






r(6 — + n q ) 
•^0 



r(6 — + ^ I ) 

Tq 0 



( 120 ) 



( 121 ) 



( 122 ) 



Here are the initial direction cosines of the 

ray at Xq , y^, and are the Cartesian coordinates 

corresponding to rQ,0Q. The parameters 6 and n are given 
by Marchand as 



n = sin 6/sin (123) 

0 = cos 9 - n cos'pg (124) 

It should be noted that Equation (117) may become singu- 
lar for certain rays where p or p^ become very close to zero 
or T . This singularity may be more easily seen in Equation 
(123) where n becomes indeterminate as both 9 and p^ approach 
zero and/or tt . In practice these conditions occur when 0^ 
is located either far out in object space, coincident with 
3, or on the image side of the lens; see Figure 8. Positioning 
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Figure 8. Excluded Regions for Center of 
Symmetry Due to Singularities 
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the center of syiranetry at these locations has therefore 
been avoided. 

B. ASSUMPTIONS 

The aforementioned assumptions for the HIN case also 
apply here. Additionally, it must be assumed, as depicted 
in Figure 7, that the GRIN lens could or soon may be fabri- 
cated from a sphere of dielectric material with the required 
spherically symmetric parabolic gradient. 

C. GRIN LENS DESIGN PARAMETERS 

In the HIN lens, the available design parameters are 
basic. These include; F, R, T, a, U, and N 2 . Parameters 
available for varying the design of the GRIN lens, however, 
include those of the HIN case but expand the index of refrac- 
tion variable N 2 into 0 , a, +b, and -b. These additional 
lens design parameters greatly expand the lens designer's 
power to bend radiant energy to his will. 

D. GRIN LENS DESIGN 

The GRIN lens design procedure, although paralleling 
that of the HIN case, is somewhat more complicated in that 
the rays are now curved and the index of refraction varies. 
Accordingly, the same design process is used but with more 
intermediate calculations required. 

The additional calculations arise since the intercept 
of the GRIN ray with the surface tangent cannot be solved 
in closed form, and an iterative solution must be used. 
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1 




The reader will note that Equation (117) will yield r if 
(Q-Qq) is known. Hence the iterative procedure is to 
"guess" (9-9 q) based on the HIN coordinates which 

are calculated as in the homogeneous lens. Thus 



0 



H 




(125) 



in Figure 9 is used in Equation (117) to find point p. The 
tangent to the ray path at point p is now extended to inter- 
cept the surface tangent again using the homogeneous inter- 
cept relations to find The prime superscripts 

indicate successive iteration values. The distance d is 
employed as a measure of the error of point p. If d is not 
within an acceptable margin then 9 ' is calculated based on 
x^H'Yxh procedure repeated to find p' in Figure 9. 

If d' is not within allowable error then the iteration con- 
tinues until it is acceptable. In practice, this iteration 
procedure has proved to be extremely rapid, rarely requiring 
more than three iterative steps before converging. Slight 
modifications, however, must be introduced to handle a posi- 
tive gradient. Furthermore, if 0^ is located outside the 
lens proper, a decreasing angular increment must be subtracted 
from to ensure that each radial vector r^^, r^^, ... inter- 

sects the ray during iteration. These modifications to the 
intercept iteration procedure are recorded in the program 
listing for program GISL (for Gradient Index Seeker Lens) . 

Refer to Appendices A, B, and C for a full description of GISL. 
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Figure 9. Grin Ray Intercept Geometry 



I 




Figure 10. Expanded View, Intercept Geometry 
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I 



with the overall iterative intercept procedure now being 



clear, the mathematical details of the GRIN lens design 
follow. Referring again to Figure 9, the GRIN ray is traced 
as in the HIN design procedure to obtain the coordinates 
X 2 (J) ,y 2 (J) on the inside surface in the meridian plane. Now, 
before Snell's law can be used, the index of refraction must 
be found at these coordinates. Using Equation (113) 



n 



2 




(125) 



where 



'J (x^(J) + OB) ^ + y 2 (J) ^ 



(126) 



and where 



R, 



Z 



y (x^(l) + OB) ^ + y 2 (l)^ 



(127) 



Furthermore, by inspection it can be easily seen that 



TT - U' - 9 



0 



(128) 



U' is identical to the HIN case^ and 



^Here Equation (125) must be substituted for N 2 in 
Ecuation ( 12) . 
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(129) 



0 



tan ^ [ 7— t rrrr 

Y2(J) + OB 




0 



In Equations (126) through (129) OB is the line segment from 
Og to B and is defined as a positive quantity to the left 
of B for algebraic ease of manipulation. By siabstitution 
into Equation (109) the scalar constant, e, of the ray can 
now be found. Furthermore, since angle U' is now known 
coordinates ^^Y te calculated using Equations (14) 

and (15) or Equations (24) through (28) . 

For the first ray defining the edge of the lens. Equa- 
tions (14) and (15) are employed, and, since no iteration 
is required, angle (9^-0^) is used in Equation (117) to find 
r immediately. Here 



X 



IH 



+ OB 



COS 9 



T 



(130) 



r, 



H 



and 



r. 



H 




(131) 



Thus 




r cos 0 



T 



OB 



(132) 



and 
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(133) 



(1) = r sin 0^ 

are the coordinates of the first point on the outside 
surface. Since e is known. Equation (119) is employed to 
find the angle ip . Angle c in Figure 10 is employed to 
translate ip into the ray direction with respect to the lens 
axis by the relation 



; = 7T - (^]; + 0^) 



(134) 



Therefore, by substituting c for U' in Equations (19), (20), 
and (21), DYDXN(l) and DYDXT(l) may be found. Here again, 
in Equation (19) must be replaced by n 2 as given by Equation 
(113) evaluated at x^^ (1) ,y^ (1) . 

Now that the first point on the outside surface is known, 
along with the surface tangent, the remainder of the K number 
of rays may be processed to yield the balance of the outside 
surface. Each siobsequent intercept between ray and tangent 
must be iterated. Thus, unlike the first ray, once ^1H' 

r^, r^, 'll, and ^ are known, ' and ' are ascertained by 

the substitution of the coordinates of point p for X 2 (J)/y 2 (>J) 
in Equations (24) through (28) with U' replaced by c. The 
coordinates of point p are 



X = r cos 0™ - OB 

p p T 



(135) 



y = r sin 0_ 
P T 



(136) 
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Therefore 









where 



A = cot(I^ + U) cot <; [y^ + Xptan;] 
B = y^(J-l) - x^(J-l) cot[Ij^ + U] 

C = 1 + cot (I^ + U) cot ? 



and 



= sin 



,-l J sin (C-U) 

( (cos(c-U) - N^/n^] ^ + sin^ (C-U) 



Note that in Equation (140) is found from Equa 
evaluated at r . 



Additionally, 



X 



IH ' 



cot Xptan^] 



Now the error, d, may be evaluated as 



d 







2 



( 137 ) 

(138 

(26) 

(139) 

(140) 
ion (113) 

(141) 

(142) 
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Here d is compared to 1 x 10 . if d is larger than this 

value then the entire procedure is repeated by substituting 

^IH '^IH ^IH'^IH forth. Once the error criteria 

are satisfied 

III* • • 

x^(J) = Xp (143) 

til* • • 

Yl(J) = Yp (144) 

and the next ray is processed. Correlation between GRIN 
(with b set equal to zero) and HIN design procedures run 
with identical parameters has shown agreement to the fifth 
and sixth decimal places. 

E. SKEW RAYS IN GRIN 

GRIN skew rays are handled analogously to the homogeneous 
case with the same coordinate transformation from grid plane 
to global coordinates being required. It is only after the 
initial directions cosines of the GRIN skew ray K', L', M' 
are found that the differences between GRIN and HIN appear. 
The only exception to this being the use of Equation (113) 
in Snell's law for refraction at the interface. Since the 
GRIN skew rays display curvature in a plane through 0^, K', 
L', and M' are constantly changing until intercept with the 
inside conical surface. Therefore, not only must the plane 
of the skew ray be analytically described, but the final 
values of K', L', 14' must be found. Due to the nature of 
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GRIN rays, the procedure for finding the ray intercept with 
the conical inside surface is different from both the HIN 
case and the iteration procedure employed in the meridian 
plane due to the multiplicity of the geometry encountered. 
The Newton- Raphson iteration routine has been found to be 
ideal for this purpose. 

To begin, the magnitude and direction of the initial 

radius vector r^ from 0 for x- / Vni must be ascertained. 

0 s 0^00 

The magnitude is given by 

1/2 

Tq = [{Xq + OB) ^ + y^ + Zq] (145) 

Therefore the unit vector in the direction of the intercept 
is 



r 



0 



Xq + OB 





(146) 



The plane of the ray may be fully described by the vector 

/N /N 

normal to the plane. Two vectors, r^ and R lie in the plane 

/N 

of the ray. R is the unit vector in the initial direction 
of the ray after refraction and described using direction 
cosines as 



R = K’i + L'j + M'k (147) 

Thus, the plane of the ray may be described by the cross 
product 
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( 148 ) 



= r. X R 

Pq 0 



or 



N =N i+N j+N k 



Ox ^Oy 



Oz 



where 



= r_ M' - r_ L' 

pqx oy Oz 



N = r K' - r* M' 
?Oy 



and 



N 



Oz 



r- L’ - r- K' 
Ox Oy 



(149) 



(150) 



(151) 



(152) 



In Ecuations (150) , (151) , and (152) r» , r„ , and r. 

Ox Uy Oz 

/N 

refer to the x, y, and z components of rQ in Equation (146) . 

/N /N 

Furthermore, the angle between r^ and R may be found from 
the dot product as 



'I'd = 

Now the scalar invariant e may be found. Substituting known 
values into Equation (109) : 
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^0 2 

= Ua + b(g-) 1 r„ sin 4,^ 

2 



( 154 ) 



where is unchanged from that found during previous calcu- 
lations for the lens shape by Equation (127) . 

With the foregoing groundwork established, the intercept 
of the ray with the inside surface may be calculated. The 
Newton-Raphson iteration scheme requires the calculation of 
the radius vector from 0^ to the cone by geometrical methods 
and the radius to the ray by GRIN theory. The difference 
between the two radii is then divided by the difference 
between the derivatives of the two functions. The resulting 

quantity is subtracted from the trial angle, 0 , in the plane 

P 

of the ray, to give a new trial angle 9 ' . The process is 

P 

continued until the difference between radii is less than 
1 X 10 The first trial angle is measured to a reference 

HIN intercept as if the material were homogeneous since the 
actual GRIN ray curves only slightly. The coordinates of 
this HIN intercept point are designated 

are derived using the HIN equations as before. To obtain 

the first trial angle, the scalar product between rg and 

r.„ is used. Here 
iH 



r. 



H 



[(x.^ + OB)^ + + 2 -h 



1/2 



(155) 



and 
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+ OB 



(156) 












i + liH 3 ^ ^ 









Thus , 



9p ' + ^Oy>^iHy + 

where ^iHz '^ ' ^ components 

of respectively. 

In Equation (117) 9 is siibstituted for (9-9») to yield 

P 0 

r as required by the iteration procedure. The geometrical 
radius, r^, is not so easily acquired. First, note that the 
equation of the plane of the ray inside the lens is given by 



Np0x<='-='0> ^ NpOy'y-yo> ^ %0z‘"-"0> = “ 

Secondly, the equation of the conical surface is given 
by Equation (75) . The combination of the plane of the ray 
and the cone yield the loci of possible intercept points on 
the inside surface. In Cartesian coordinates, the sum of 
Equations (75) and (158) is 



x(N -X tan oc] + y[N +y] + z[N. +z] ~N X- 
pox ^ poy ^ pOz pOx 0 



-N„v. -N^z» = 0 

pOy^O pOz 0 



(159) 



Equation (159) must be transformed into coordinates r and 0 
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in the plane of the ray. The transformation is made possible 
using Equations (120) through (124) . Upon substitution, 



and after solving for r^, Equation (159) becomes 



— 1 - \ —1 + J - 

2a ~ V 2 a 

2 4A2 2 



where 



2 2 2 2 
^2 " ®1 ^1 ■ “ 



= B,N + C,N + A,N + 2A, OB tan a 
2 1 poy 1 poz 1 pox 1 



■^2 “ °®^pox 



+ N Z- 
poz 0 



and 



A, 



— (x + OB) + nK' 
^0 ^ 



B. 



oy. 



+ nL' 



oz, 



+ nW 



Here, 6 and n are found from 



n = sin 6p/sin 



(160) 

(161) 

(162) 

(163) 

(164) 

(165) 

(166) 

(167) 
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and 



5 = cos 0 - cos 

P 0 



( 168 ) 



which follow from Equations (123) and (124). The plus sign 
in Equation (160) yields the correct values. Now that r of 
the rav and r of the surface are known, the derivatives of 

g 

r and r with respect to 8 at 9 must be foxind. It can be 
g P 

shown that for the ray 



dr 

de 



— cos [-2 £ 9 + sin (A,) ] 

- P 3 

2e 



(169) 



where 



„ 2 , 2 

^a^ + 4be^/R2 



(170) 



The derivative of r is somewhat more complicated. With 

g 

persistence, however, it can be shown that 



dr 

d9 



dB. 



dA„ , bI C- -1/2 
^^^^2 "dl ®2 "dF^ ^ ^^4^ 



(171) 



where A^, B^, and C 2 are given by Equations (161) through 
(168) and 



1 2 ‘^^2 



2 *^2 
“2 2 d0 ’ 



^2 *^2 



2 d9 



(172) 
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Therefore, the r and 9^ to the intercept are found by 
iteration of revised trial values 



(r - r ) 



9 



9 



2 . 



( 173 ) 



P 



P 



(— 



d9 d9 



In practice the quotient of differences in Equation (173) 
is reduced by a factor of 1.3 to slow convergence and provide 
stability. The number of iterations, however, rarely exceeds 
five . 

Values for r and 9 are now transformed into Cartesian 



coordinates x^, y^, by substitution into Equations (120) 
through (124) . Next, the values of the direction cosines 
K', L', and .M' at the intercept are needed. The angle 4> 
between the radius vector and the tangent to the ray at 
intercept may first be deduced from the scalar invariant, 
e. Hence 



Ifhere e is known, is evaluated at rp^ is the 

radius to the intercept as found by the iteration above. 
Three constraints on the direction cosines may be written. 
These are: 

1) The scalar product of the radius vector, ^p^, and the 



instantaneous ray direction vector, R, at intercept. 

2) The scalar product of the normal to the plane of the 

/N /> 

ray N ^ and R. 



P 




(174) 
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3) The sum of the squares of the direction cosines must 



sum to unity. 

Mathematically/ the above constraints are written as 
(x. + OB) y. z. 

— K' + L' + M' = cos ijj (175) 

PI PI ^PI 



^POx^' 



N . L' 
pOy 



N ^ M' 
pOz 



= 0 



(176) 



2 2 2 
K'^ + L'^ + = 1 



(177) 



To find K', L', and M', Equations (175), (176) and (177) are 
solved simultaneously. It can be shown that the solution 
leads to 




(178) 



L' = 



D 



B„M' 

D 



(179) 



K' 



= (1 - 



2 2 
L' - M'^) 



1/2 



(ISO) 



where 



"pOx 



No o 

+ (^)2 ^ 1 
'’pOx 



(181) 



94 



(182) 



B 



6 




C 



6 



U + (^E2£)2]a2 - 1 



pOx 



(183) 



and 



A, 






(184) 



5 





B 



5 




(185) 



Furthermore, in Equations (184) and (185), , B^ , and 

are the coefficients of K', L’ and M' in Equation (175) . 

Now that the direction cosines of the ray are known at 
the point of intersection with the inside surface, the index 
of refraction, n^ , is computed by substituting rp^ into 
Equation (113) . The remainder of the skew ray trace to the 
image plane is identical to the homogeneous procedure. 
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IV. LENS PERFORMANCE PARAMETERS 



The function of the seeker lens is to focus electromag- 
netic energy either reflected from or emitted by the target 
onto a detector. Angular displacement of the target with 
respect to the missile body axes as well as target angular 
rate information are both desired outputs from the seeker. 
Hence, it is not only important just to be able to detect 
the target by focusing energy into a spot on the detector, 
it is equally important that this spot be as small as possi- 
ble to enable the precise position of the spot on the detec- 
tor to be discerned. 

The ability of a lens to focus an object to a small spot 
does not guarantee the quality of the image. For a FLIR 
optical system extensive effort is expended to obtain an 
image with minimum aberration [10]. Seeker optics, however, 
are generally non- imaging devices where the pressure of the 
different aberrations does not detract from the function of 
the seeker as long as a tight image is maintained [11] . 

Accordingly, the most important parameter by which seeker 
lens performance is judged is that of spot size at different 
obliquities. Since the image found at the focal point is not 
necessarily circular nor equally dense, the standard devia- 
tion in the y and z directions with respect to image centroid 
is used to define spot size. 

Therefore, once the coordinate pairs of all the skew rays 
have been calculated in the image plane, the first step in 
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the analysis of lens performance is to calculate the image 
centroid [12], If the spot diagram is composed of N rays. 



Since there is symmetry about the x-y plane, the z coor- 
dinate of the centroid will always be zero. 

The next step in finding the spot size is to find the 
standard deviations of the spot diagram in the y and z 
directions . This is accomplished by summing the squares of 
the differences of the intercept coordinates with respect 
to the centroid and then dividing by the number of rays. 

The standard deviations then, are given by 



the centroid location, y is 

c 



V 

' c 




(186) 



a 



z 



2 




(187) 



and 



a 




(188) 



The spot size, o^, is now defined by 



a 



r 




(189) 
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Of further interest in appraising lens performance is 
the energy density of the image as a function of radius from 
the centroid. In nondimens ional form, this is simply the 
number of rays in the spot diagram within a succession of 
circles of increasing size overlaid about the centroid. Here, 
each ray is assumed to carry a unit amount of radiant energy. 
This type of plot facilitates the comparison of different 
lens designs by detailing the distribution of energy within 
each image. Clearly, it is desirable to have as much energy 
as possible concentrated very close to the image centroid. 
Between two lenses with equal spot sizes, the preferred lens 
has more energy concentrated within a smaller radius. 

Every ray, however, does not deliver an equal amount of 
energy to the focal plane. It is prudent, therefore, to 
include as a performance parameter the average ray intensity. 
Again, for N rays, the intensity I of each ray is summed and 
normalized by N to yield 



I 

a V 



I I 
■=1 



( 190 ) 



Finally, each lens design is checked for "failed rays". 
The reader will recall that these rays fail to intercept the 
image plane due to total internal reflection, total external 
reflection, or failure to intercept the inside surface within 
the bounds of the lens. Hence, a lens design with fewer 
"failed rays" or no "failed rays" at all is a preferred lens. 
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Although there are many other performance criteria by 
which lenses are compared, the foregoing parameters are more 
than sufficient to judge the merit of preliminary seeker 
lens designs. It should be noted, however, that notwithstand 
ing the fact that the image centroid and standard deviations 
were used as stepping stones to obtain image spot size, they 
have significant meaning of their own. The standard devia- 
tions and inform the lens designer as to the horizontal 
and vertical spread of the image. Image centroid location, 
y^, at increasing obliquities is of obvious importance since 
excessive displacement will cause the image to miss the 
detector entirely and would dictate the necessity for a 
second lens element to dampen the movement. Furthermore, 
since Line of Sight (LOS) measurement accuracy to the target 
is highly dependent upon the linearity of y^ as a function of 
the lens tilt angle, lens designs which exhibit a greater 
degree of such linear behavior are the preferred designs. 

To summarize, it is sufficient to note that although 
spot size is the most important of the performance parameters 
every other parameter has a significant impact on the 
performance of a particular lens. 
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V. RESULTS FOR THE HOMOGENEOUS LENS 



The performance of the homogeneous lens is presented 
primarily as a comparison with which to compare the performance 
of the GRIN lens. Here, the relationship of spot size to 
increasing lens obliquity and lens thickness are presented 
as well as the image centroid movement as a function of 
obliquity. Additionally, the reader is introduced to the 
four basic computer plots used to display the results: lens 

shape, object plane with superimposed skew ray grid, image 
plane or spot diagram, and energy density. These plots were 
generated on the VERSATEC Plotter using arrays of data points 
produced by program GISL on the IBM 3033 mainframe computer. 

To begin. Figure 11 shows the homogeneous lens shape. 

The first of the four basic plots, the lens shape plot, pre- 
sents the lens side view in the meridian plane. The outer 
surface (curved) and the inner surface (cone) are constructed 
by connecting the points x^(J) , y^(H) and X 2 (J) ,y 2 (J) / 
respectively, by straight lines. At the apex of the lens 
about the lens axis is a trapezoidal region which represents 
the opaque nose area. In the legend are listed the lens 
design parameters and the significant calculated dimensions 
of the lens. Since all linear dimensions are implicitly 
normalized with respect to the maximum inside radius, R, 
the lens may be scaled up or down by multiplying each dimen- 
sion by a factor of Beginning at the top, the 

parameters listed in the legend are explained in Table 1. 
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LENS SHAPE 




F/« = F/2 
F = 4.000 
OSYMB = 1.000 
R = 1.000 
ALPHA = 0.7853975 
T = 0.050 
U = 0.00000 
I = 1008 
N1 = 1.0000 
N3 = 1.0000 
A = 2.250000 
B = 0.000000 
BEFF =0.0000000 
7. INDEX = 0.000 
7. INSIOE = 0.0000 
7. OUTSIDE = 0.0000 
7. ACROSS =0.0000 
STATION A =-0.4417198 
GAMMA = 0.7853968 
TSURFL = 1.7415296 



x-lems axis 



Figur6 11- Homogeneous Lens Shape for N 2 1.5 
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TABLE 1 



Explanation of Lens Shape Plot Legend 



PARAMETER 

F/# 

F 

OSYMB 

R 

ALPHA 

T 

U 



N, 



N. 



B 

®EFF 
% Index 

% Inside 
% Outside 



TYPE OF PARAMETER 
DESIGN 
DESIGN 
DES IGN 



DESIGN 

DESIGN 

DESIGN 

DESIGN 

DESIGN 

DESIGN 

DESIGN 

DESIGN 



DESIGN 

CALCULATED 

DESIGN 

CALCULATED 

CALCULATED 



MEANING 



F nuit±>er. F/# = F/2R 

ic 

Focal length from B . 

OB. Line segment O to B 
(positive to left) 
Immaterial in HIN 



Maximum radius of cone 



a — cone half-angle, radians 
★ 

Edge thickness 

Incident ray offset angle 
(design) *, radians 

Number of iterations. 

I+l = number of rays 

Free stream index of 
refraction 

Index of refraction of 
interior lens cavity 



a in 



n^(r) = '/at b (r^/R^ , 



gradient refractive 
index function, 
in HIN. 



Nj = /a 



b in n 2 (r) . Zero in HIN 

2 

b effective = b/R^ 

Percent change in n_(r) 
from r = 0 to r = R^ 

Percent change in n„(r) 
along inside surface from 
lens axis to edge 

Percent change in n 2 (r) 
along outside surface 
from opaque region to edge 
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Table 1 (Continued) 



PARAMETER 


TYPE OF PARAMETER 


MEANING 


% Across 


CALCULATED 


Percent change in n2(r) 
across lens from lens 
axis to outside surface 
at the thickest point 


STATION A 


CALCULATED 


x-coordinate at nose of 
lens* 


GAMMA 


CALCULATED 


y — nose half-angle of 
opaque nose region, 
radians 


TSURFL 


CALCULATED 


Total outside surface 
length from Station A 
to the edge 



X 

Refer to Figure 2 for clarification 
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HIN Isns shape has a convex outer surface with maxi— 
mum thickness on axis of almost ten times the edge thickness. 
Although the lens has a good aerodynamic shape resembling an 
ogive, the outer surface is not a circular arc nor can a 
single analytical function be fitted to the array of points 
describing the surface. Note that the nose half-angle, y, 
is almost identical to the cone half-angle, a. 

All lenses have been designed with a cone half-angle of 
45° which is approximately the maximum angle for which aero- 
dynamically efficient lens shapes may be designed, consider- 
ing a free stream Mach Number not to exceed three. Without 
exception, overall lens performance is more severely degraded 
as angle a is reduced. 

Table 2 explains the legend of Figure 12, which is the 

second basic plot. Here the lens is depicted as seen from 

the skew ray grid plane. The lens tilt angle, a , causes 

P 

the equally spaced (in J) circles descriging the surface of the 
lens to appear as ellipses. In Figure 12, the grid spacing has 
been reduced from 0.1, which is normally used, to 0.3 to 
allow identification of individual rays for correlation with 
the image plane spot diagram. Although the small number of 
rays used is not sufficient to given an accurate definition 
of spot size, the number is sufficient to describe where 
rays in the object plane are being focused in the image plane 
by the HIN lens. The skew rays in Figure 12 have been numbered 
in the order in which they were processed. Actually, only 
rays 1 through 19 were actually traced; 20 through 33 are 
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TABLE 2 



Explanation of OBJECT PLANE Plot Legend 



PARAMETER TYPE OF PARAMETER 



MEANING 



ALFAP 



Analysis 




Rays 



Analysis 



Total number of rays 
processed 



Failed Rays 



Analysis 



Total number of rays 
failing to pass through 
the lens* 



Failed Rays are indicated on the plot by a diamond 
superimposed on the grid location of the ray. 

mirror image skew rays. Ray 30 corresponds to ray 14, for 
example. Here, the lens has been tilted by 0.4 radians, or 
22.9 degrees, and 33 rays have been processed of which none 
have failed to intercept the image plane. The staircase 
pattern has been added in this case in order to show the 
resulting distortion present in the image plane (Figure 13) . 

The Spot Diagram in Figure 13 is an example of the third 
basic computer plot; see Table 3. Unlike most Spot Diagrams, 
this example has the individual rays numbered for comparison with 
Figure 12; also the resulting distorted staircase pattern 
is sketched. By cross-referencing individual rays between 
Figures 12 and 13, it is possible to recognize where certain 
areas of the lens are focusing rays in the image plane. Rays 
1, 4, 20, and 22 about the opaque nose region form a coma 
tail which contributes most of the image spread. Rays 33, 
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Figure 12 . Example Object Plane 
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Figure 13. Spot Diagran for HIN Lens Design Shown 
in Figure 11 
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TABLE 3 





Explanation of Spot 


Diagram Leqend 


PARAiMETER 


TYPE OF PARAMETER 


MEANING 


ALPAP 


Analysis 


Up, lens tilt angle, 
radians 


SPOT SIZE 


Performance 


a^, see Equation (189) 


Y-CENTROID 


Performance 


Yet y-coordinate of 
image centroid; see 
Equation (186) 


SIGMA Y 


Performance 


Oy, y-standard deviation; 
see Equation (188) 


SIGMA Z 


Performance 


02 f z-standard deviation; 
see Equation (187) 


RAY COUNT 


Performance 


Number of rays striking 



the image plane 

18, and 19 from the bottom portion of the homogeneous lens 
are imaged at the top and are widely separated from the core 
of the image. In general, the regions of the lens which 
have been found to contribute the bulk of the widely spaced 
rays are the immediate nose region and the lower portion of 
the lens . 

Furthermore, the distortion present in the image of the 
staircase pattern clearly shows that regions closest to the 
nose yield the greatest distortion. Horizontal lines are 
switched end for end and tilted approximately 45 degrees. 

The upper portion of the lens performs the best. Rays 
3, 7, 11, and 15 in the meridian plane are focused on the 
y^^ axis in a fairly tight region; rays 10, 14, 17, 27, 30, 



and 32 about the upper periphery in the object plane are all 
focused within the image core. 

Table 3 explains the legend of the spot diagram. It is 
seen that the standard deviation in the y— direction is 
approximately five times that of the z-direction. This 
elongation of the image is not readily evident in the spot 
diagram since the ordinate and abscissa have not been plotted 
with equal increments. This results from the great disparity 
between Spot Diagrams of the various lens designs studied. 

It is important, therefore, for the reader to take careful 
note of the relative sizes of the y^,^ and z^,, axes. 

Of primary importance is the spot size in Figure 13. The 
value for spot size is adversely affected by the poorest 
performing regions of the lens. Were it not for these 
errant rays, the spot size would be considerably smaller. 

The image intensity pattern is benefitted, however, by the 
fact that the rays spread the farthest from the centroid 
contribute significantly less energy per ray than those 
being focused in the core of the image. Table 4 lists the 
relative intensities of the primary skew rays plotted in 
Figure 13. It is seen that the high intercept angles experi- 
enced by the rays closest to the bottom of the lens and, to 
a lesser extent, those near the nose, result in higher 
reflectivity and lower transmission through the lens. 

In the legend of Figure 14, the average of ray intensity 
is given. The Spot Diagram Energy Density distribution may 
be seen at a glance. The fraction of energy (number of rays) 
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TABLE 4 



Skew Ray Intensities of HIN Lens 



a) Numberical Order 



Ray 


Intensity 


1 


0 . 746 


2 


0 . 789 


3 


0.921 


4 


0 . 89 2 


5 


0.877 


6 


0.811 


7 


0.920 


8 


0.917 


9 


0.904 


10 


0 . 870 


In order 


of Descending 


Intensity Ray(s) 


0.921 


3 


0 .920 


7 


0 .919 


11 


0.918 


12,15 


0.917 


8 


0 .916 


16 


0.910 


13,17 


0.904 


9 


0.89 2 


4 



Ray 

11 

12 

13 

14 

15 

16 

17 

18 
19 



Intensity 



Intensity 


Ray (i 


0 . 887 


14 


0.877 


5 


0 . 870 


10 


0.811 


6 


0.789 


2 


0.746 


1 


0.640 


18 


0.576 


19 



Intensity 

0.919 
0.918 
0.910 
0 . 887 
0.918 
0.916 
0.910 
0.640 
0.576 
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Figure 14 . Nondimensional Encircled Energy Plot for 
HIN Lens Design of Figure 11 
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shown as a function of radius from the Spot Diagram centroid 
is normalized with respect to the total number of rays 
traces (Rays, in Figure 12) regardless of whether all of 
the rays successfully intercepted the image plane. Thus 
the fourth basic plot may be used in conjunction with the 



Spot Diagram to further define the image concentration with 
respect to the centroid. 

The response of the HIN lens to increasing tilt angles 
is given in Figures D-1 through D-23 where the index of refrac- 
tion of the lens has been set at 1.5. Whereas the spot size 

at a = 0.0 is very small, that at a =0.7 radians is quite 
P P 

large at 67% of the lens radius. Figures 15 and 16 summarize 

spot size and centroid locations for the lens. Spot size 

growth, although somewhat irregular at the higher angles, 

is pronounced. Furthermore, the centroid movement is seen 

to be an approximately linear function of a until 

P 

a = 0.3 and easily exceeds the radial dimension of the lens 
at higher tilt angles. 

Figure 17 shows that the HIN lens may be slightly improved 
by increasing the edge thickness, T, by a small amount. The 
lens may be otherwise tuned to improve performance at certain 
tilt angles by designing the lens with U slightly greater 
than zero. These performance improvements are practically 
insignificant, however, and neither lens tuning by the param- 
eter U nor T produce improvement across the spectrum of tilt 

angles. Instead, an improvement at one a usually has resulted 

P 

in a degradation at others. 
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SPOT SIZE VS. LENS TILT ANGLE 




(sniQva swa) 321s loas 



Figure 15. HIN Lens Spot Size versus a 
for = 1.5 P 
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MOVEMENT OF SPOT CENTROID 




Figure 16 . HIN Lens Centroid Movement versus 

a at N_ = 1.5 
P 2 
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ALFAP (RADIANS) 



SPOT SIZE VS. EDGE THICKNESS 




(sniGva swa) 2zis loas 



Figure 17. HIN Lens Spot Size versus Edge Thickness 

at a =0.3 Radians and N_ = 1.5 
P 2 
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EDGE THICKNESS 



Finally, the homogeneous lens is shown with index of 
refraction of three at the intermediate tilt angle of 0.3 
radians in Figures D-24 through D-27. The lens shape required 
to accommodate the higher index of refraction is seen to be 
thinner and displays less outside surface curvature than 
the HIN lens with N 2 = 1.5. Spot size is significantly re- 
duced. Since infrared lenses such as germanium generally 
have fairly high refractive indices, improved lens per- 
formance at these higher values is encouraging. 
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VI. GRIN LENS RESULTS 



The performance of the GRIN lens is similar in many 
respects to the homogeneous lens. The relationship between 
characteristic regions of the GRIN lens, such as the lower 
lens portion, and where these regions image bundles of rays 
is identical to the HIN lens as depicted in Figures 11 and 
12. The growth of spot size with increasing and the 
respective movement of image centroid typical of the homogene 
ous lens is clearly displayed by the GRIN lens as well. 

The measure of the superiority of the GRIN lens, there- 
fore, lies in the successful correction or improvement of 
the deficiencies seen in the HIN lens. Here, the reduction 
of spot size is of primary concern. 

As a modest example of the ability of the gradient refrac 

tive index to reduce spot size, the GRIN lens design shown 

in Figure 18 is examined. This lens is very similar in shape 

to the homogeneous lens with N 2 = 1.5. Note, however, that 

unlike the HIN lens, the nose half-angle, y, is slightly 

larger than that of the cone, indicating more outside surface 

curvature. The object plane for this lens, at = 0.4, is 

shown in Figure 19. Hers, the error in focusing rays present 

in the image plane has been superimposed over the grid plane. 

y.... coordinate error contours (with respect to the centroid) 
IM 

are shown in the left half of the plane and z^^^ contours in 
the right. These contours vividly show that the largest 
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LENS SHAPE 




F/« = F/2 
F = 4.000 
OSYMB = 1.000 
R = 1.000 
flL^HR = 0.7853975 
I = 0.050 
U = 0.00000 
I = 1008 
N 1 = 1.0000 

43 = 1.0000 
fl = 2.950000 
B = -0.427500 
BEFF =-0.0855000 
7. INOEX = -10.000 
7. ^NSIOE = -3.2397 

7. OUTSIOE = -9.2302 

7. ACROSS =1.0913 
station a = -0.43525^^ 
GAMMA = 0.8068293 
^SURFL = 1 . 7399152 



X-LENS QX:S 



Figure 18. Example GRIN Lens Design with 10% 
Negative Gradient at OB = 1.0 
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Figure 19 . Example GRIN Lens Object Plane. Contours in the Left Half Depict 

the Difference Between the y Coordinate in the Image Plane and Image 
Centroid. The Right Hand Contour Shows the Corresponding z Displace 
ment in the Image Plane. 



errors in both y and z stem from the nose and bottom of the 
lens. The best performance is contributed by the upper 
central region. The reader will note that negative errors 
in z occur below the y^ = 0 plane of the lens and that nega- 
tive errors in y straddle the = 0 plane. 

The corresponding intensity contours are displayed in 
Figure 20. Of significance here is that regions which perform 
relatively well in y and z error also perform well in trans- 
mitting energy. Regions of the lens which are characterized 
by relatively high angles of incidence, therefore, perform 
the poorest and have the most to gain from better combinations 
of gradient index and Og 

The Spot Diagram, Figure 21, displays an image pattern 
virtually identical with the HIN lens in Figure 11. The 
GRIN spot size, a^, however, is slightly smaller by approxi- 
mately 2% and both and are correspondingly smaller 
which indicates superior performance by the GRIN version. 

Image centroid for the GRIN is displaced further than that 
for HIN. The average intensity transmitted by the GRIN lens 
compares favorably at 0.88 as opposed to the average inten- 
sity of the HIN at 0.85. Further evidence of superiority 
of the GRIN lens is seen in the form of a steeper slope to 
the Encircled Energy Plot (Figure 22) than that of the HIN 
lens in Figure 13. 

Although the example GRIN design displays marginal 
superiority, some GRIN designs do not. In particular, it 
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Figure 20 . Intensity Contours of the Example GRIN Lens Shown in 
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Figure 21. Example GRIN Spot Diagram at Up = 0.4 
Radians for GRIN Lens Design Shown in 
Figure 18 
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SPOT DIAGRAM ENERGY DENSITY 
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Figure 22. Example GRIN Encircled Energy Plot 
Corresponding to Figure 21 
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has been found that as the center of symmetry of the index 
function, 0^, is moved further away from the interior of the 
GRIN lens, performance deteriorates. 

Two series of GRIN lens designs have been examined at 
the intermediate tilt angle of 0.3 radians in order to define 
the spectrum of lens performance. Both design series explore 
the GRIN lens for 0^ inside, outside, and far outside the 
interior of the lens. The first series. Figures E-1 through 
•C'~92 , are GRIN designs in the "low range" of refractive index 
with the index parameter "a" set at 2.25. This first set is 
compared to the HIN lens with N^ = /a = 1.5. The second, 
or "high range" series is for a = 9.0, as compared to the 
HIN at N 2 = 3.0, and are displayed in Figures F-1 to F-128. 
Both the high and low range lenses exploit gradient changes 
of 5, 10, 25, and 50 percent both positive and negative, 
where possible. A negative 50% gradient in the low range 
is, of course, not possible since an index of below 1.0 
would result. The center of symmetry was not located at 
X ^ 0 due to the aforementioned singularity encountered at 
very small angles. 

The lens shapes of the GRIN designs shown vary widely. 

The outside surface may be either convex or concave and ex- 
hibits a thinner profile at higher values of refractive 
index. Although all of the resultant GRIN lens shapes are 
superior aerodynamically to the hemispherical seeker lenses 
currently in use, the convex version of the lens has more 
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obvious applications for a Sidewinder type tactical missile, 
whereas the concave lens shape is more applicable to the 
diffuser of a ramjet with nose inlet. 

Object plane diagrams show that for lens designs exhibit- 
ing smaller spot size, fewer or none of the skew rays failed. 
This fact underscores the success of GRIN in controlling image 
deterioration contributed by the nose region of the lenses. 
Although the lower lens and nose regions still create the 
largest of the image spread, the extent is reduced. 

Figures 23 and 24 are contour diagrams summarizing 
spot size performance of GRIN lens designs in the low and 
high range respectively. From these diagrams, it can be seen 
immediately that the best performance is obtained from lenses 
having the center of symmetry inside the lens. From this 
fact, it may be deduced that large changes in refractive 
index are desired along the surface of the lens rather than 
across. Except for isolated regions at -5%, increasing 
gradients produce smaller spot sizes. The positive gradient 
in the high range, however, has an almost constant spot size 
from +5% to +50% where an improvement of only 0.6% is seen. 

In both the low and high ranges the positive gradient at 
50% and at OB = 0.05 has proved to be the best performing 
combination. Again, as in the HIN lens the high range 
exhibits the smaller spot size. 

The best GRIN lens is now further examined over the range 
of obliquities up to 0.8 radians. Furthermore, since it has 
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Figure 23. Contour Plot Summary of GRIN Lens Spot Size Performance 
for a =0.3 radians, a = 2.25 
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Figure 24 . Contour Plot Summarizing GRIN Lens Spot Size Performance 
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been determined that placing the focal point at x = 2.0, 

improves spot size performance at lower values of a , another 

P 

series of plots is given for F/1. With the focal point at 
X = 1-5, severe degradation results; for F = 3.0 only modest 
differences with F/2 are evident. Appendices G and H show 
the performance of the "best" GRIN lens at F/1 and F/2. 

Here it may be fully recognized that the resulting be- 
havior of a "good" gradient refractive index seeker lens is 
far superior to that of the HIN lens for a <_ 0 . 6 radians; 
see Table 5. The reader will note that there are no failed 
rays even at 0.8 radians (45.8 degrees). Although centroid 
movement at F/2 is slightly greater than the HIN lens, by 
changing to F/1, this can be corrected; see Figure 25. The 
F/1 version, moreover, displays a range of spot sizes signi- 
ficantly smaller than that for the F/2 lens below a =0.54 

P 

radians or the HIN lens below 0.68 radians; see Figure 26. 

That the lens displays more desirable performance with a 
shorter focal length is a surprising but highly desired re- 
sult. Since seeker optics systems are volumn limited, any 
reduction in the lens focal length is beneficial to the final 
design and packaging requirements. 

It is significant for the reader to note that centroid 
movement beyond a value of ±1.0 exceeds the physical dimen- 
sions of the lens radius. If the interior of the lens mount 
is of similar dimension, it follows that the requirement for 
y„,,„_ to be less than rl.O restricts the F/1 lens to approxi- 
mately 0.65 radians of tilt whereas the F/2 lens is further 
restricted to less than 0.3 radians. The superior performance 
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Figure 25. "Best" GRIN Lens Centroid Moment for 
F/Numbers of F/1 and F/2 
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Figure 26. "Best" GRIN Lens Spot Size Performance 
for F/Numbers of F/1 and F/2 
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TABLE 5 



Comparison of HIN and GRIN Lens Designs 



Performance Parameter or Feature 
Shape of Outer Surface 

Fraction of Failed Rays 

at a =0.7 radians 
P 

Linearity of Centroid 
Motion; a_ for Ay/y 
deviation^of 10% 



HIN 



Convex 



0.14 3 



GRIN 



Concave 



0.0 



0.4 radians 0.8 radians 



Spot Size for a 



a 


ZZ 


0.1 


radians 


0 . 190 


0.029 


p 


a 

r\ 




0.2 


radians 


0.247 


0.063 


y 

OL 


— 


0.3 


radians 


0.442 


0.112 


p 


a 

r> 


= 


0.4 


radians 


0.466 


0.201 


a 

n 


= 


0.5 


radians 


0.637 


0.313 


y 

OL 




0.6 


radians 


0.632 


0.490 



a =0.7 radians 


0.667 




0.688 


ir 

a =0.8 radians 
P 


N/A 




0.947 


Spot radius for 80% 
of energy § a =0.4 
radians ^ 


Does not 
attain 80% 
within 0.50 


0.28 


Average relative intensity 
of skew rays 5 = 0.4 

radians 


0 . 82 




0.44 



*"Best" GRIN lens with 50% positive gradient, 0^ - 0.05, 
in the F/1 configuration. 
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of the F/2 lens above 0.54 radians iS/ therefore/ unusable. 
Attempts to improve upon this performance by increasing the 
edge thickness, T, or by slightly adjusting the parameter U, 
as with the homogeneous lens only produced uniformly degraded 
performance in every respect. 

Although the "best" GRIN lens resulted from a positive 
50% gradient change in refractive index, it has been noted 
that this configuration was only slightly better than the 
same lens with a positive 5% change. Furthermore, since 
refractive index gradients of five percent or better have 
already been produced, it is entirely feasible that if the 
precise parabolic change could be controlled, this lens could 
be produced today. 

Despite the obvious success of GRIN in controlling spot 
size growth and image centroid movement, a penalty in the 
form of reduced ray intensity has been paid. That increasing 
spot size performance is tempered by a loss of intensity may 
be seen by comparing "good" GRIN encircled energy plots with 
that of the HIN lens. Note that the HIN lens with N 2 = 3.0 
also loses - 50% to intensity; see Figure 44. This loss in 
intensity is partially offset, however, by the GRIN lens 
with an increased number of rays transmitted and by the re- 
duced spot size. 

One drawback to the use of the GRIN lens as a self 
sufficient, single element lens is still that of the relative 
sizes of image and detector. Hence, even though GRIN has 
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significantly reduced spot size below a =0.6, that size is 

P 

still significant at tilt angles above 0.3 radians. In order 
to use the lens without a secondary focusing element or 
mirror arrangement, a large composite sensor array would be 
required. 
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VII. CONCLUSION 



Gradient refractive index materials may be employed to 
design a pointed seeker lens which exhibits optical per- 
formance far superior to that obtainable with conventional 
homogeneous optical material. A fifty percent, positive, 
parabolic gradient index with center of symmetry interior 
to the lens was found to yield the best performance although 
a five percent version of the same lens was only very slightly 
inferior; this lens may possibly be fabricated today. 

For objects off-axis by more than 0.3 radians (-17.2 
degrees) a secondary lens element may be required unless a 
large scale multiple element sensor array is employed. With 
such arrays, objects off-axis by more than 0.65 radians (37.2 
degrees) may require Cassegrainian or other mirror elements 
to compensate image movement. 
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VIII . RECOMMENDATIONS FOR FUTURE WORK 



This thesis investigated a spherically symmetric GRIN 
seeker lens with inside conical surface. Future studies 
should investigate : 

a) The elimination of the singularity problem by 
modification of the theoretical equations for GRIN 
using numerical techniques. See Reference 13. 

b) The optimization of the GRIN lens with both inside 
and outside curvature. 

c) The effect on lens performance due to an attached 
shock wave. 

d) The effect on lens performance when the object is 
no longer in the far field and wave front curvature 
must be taken into account. 

e) The performance of the lens using wavelengths corres- 
ponding to atmospheric windows. 

f) The effect of a radially symmetric gradient index 
on lens performance. 

g) The feasibility of adding an anamorphic gradient- 
index lens located on the missile body to increase 
off-axis tracking/acquisition capability to 90° and 
beyond. See Figure 1 of Reference 14. 
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APPENDIX A 



COMPUTER LOGIC FLOW DIAGRAM FOR PROGRAM GISL 

This appendix together with Appendices B and C describe 
the FORTRAN program GISL. GISL may be used to design either 
a homogeneous (B = 0) , or a GRIN lens to the user's specifi- 
cations by changing the design parameter where indicated by 
comments in the input section of the program. Additional 
software required to run GISL and plot the results are not 
described since these system-dependent procedures do not 
apply elsewhere. 

The flow diagrams which follow provide the reader with 
the information necessary to follow the fundamental logic 
of the main program and subroutines of GISL. Not shown is 
Subroutine DIRECT which has been derived from Subroutine 
DIRECT as described by Amichai [5] with only slight modifi- 
cations . 
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F.ere, the priire indicates Y-ooordinate as transfonied into the grid plane. 
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APPENDIX B 



LISTING OF FORTRAN PROGRAM GISL 
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IMAGE PLANE SPOT CIAGRAM ANALYSIS: 
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TO OBTAIN PLOTS FROM TH£ PLOTTER OP THE PRINTER t 
ISSL6 THE FOLLOWING COPLAND: 

FOR PRINTER PLOTS ENTER ‘’CHARTS PRINTER” 

FOR PLOTTER GR/PHS ENTER ’’OHARTS PLOTTER" 

END OF PROGRAM 
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APPENDIX D 



HIN LENS PERFORMANCE PLOTS 
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Shown in Figure 11 
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. Spot Diagram Corresponding to Object 
Plane of Figure D-1 
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LENS FRONT VIEW 
OBJECT PLANE 
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Figure D-3. Object Plane at a„ = 0.1 Radians for 
HIN Lens Shown in^Figure 11 
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Figure D-5. Encircled Energy Plot for the Spot 
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Figure D-6. HIN Object Plane at a = 0.2 Radians 
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Figure D-9 . HIN Object Plane at a =0.3 Radians, 
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Figure D-12. HIN Object Plane at a =0.4 Radians 
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Figure D-22. Spot Diagram of HIN Lens at a = 0.7 
Radians, = 1*5 
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D-23. Encircled Energy Plot at - 0 . 7 
Radians, N 2 = 1-5 



760258 

5 



204 



-MERIDinN PLANE 



LENS SHAPE 




F/» = F/2 
F = 4.000 
0SYM8 = 1.000 
R = 1.000 
ALPHA = 0. 7853975 
T = 0.050 
U f 0.00000 
I = 1008 
N1 = 1.0000 
N3 = 1.0000 
A = 9.000000 
8 = 0.000000 
BEFF =0.0000000 
7. INOEX = 0.000 
/. INSIDE = 0.0000 
7 OUTSIDE = 0.0000 

7. ACROSS =0.0000 
STATION A = -0. 1727498 
GAMMA = 0. 7853975 
TSURFL = 1.5313670 



X-LENS RXIS 



Figure D-24 . HIN Lens Design for 1^2 “ 
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Figure D-25. Object Plane of HIN Lens in Figure 
at a =0.3 Radians 
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Figure D-27. Encircled Energy Plot for Spot Diagram 
of Figure D-26 
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Figure E-1. GRIN Lens Shape at +50%, OB = 0.05, 
a = 2.25 
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Figure B-2. Grid Plane at a =0.3 for Lens of Figure 
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Figure E-4 . Encircled Energy of Figure E-3 
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Figure. E-6. Grid Plane at a =0.3 for Lens of Figure 
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215 



h NLhu 



SPOT DlflGRRM ENERGY DENSITY 



. . 'j 



ij . 6 



— I]. S 

j 



o.q 



j 






jm 

f ‘ 




u 



J. 2 



i 



/ 



■ I 

0.2 C . 1 



RLFRP - 0.300000 
RVERRGE Rh" INTENSI T T - 0 



,i,l. 



I I ' 



Q.2 0.3 O.q 0, 

'PGIU5 FROM IMRGE CEN'ROID 
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Figure E-9. GRIN Lens Shape at -2%, OB = 0.20, 
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Figure E-10. Grid Plane at a =0.3 for Lens of Figure 
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Figure E-12. Encircled Energy of Figure E-11 
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Figure E -14 . Grid Plane at a = 0.3 for Lens of Figure E -13 
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Figure E-16. Encircled Energy of Figure E-15 
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Figure E-18. Grid Plane at a =0.3 for Lens of Figure E-17 
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Figure E-19 , Spot Diagram for Grid of Figure E-18 
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Figure E-20. Encircled Energy of Figure E-19 
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Figure E-22. Grid Plane at a =0.3 for Lens of Figure E-21 
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Figure E-24. Encircled Energy of Figure E-23 
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Figure E-25. GRIN Lens Shape at +25%, OB = 0.20, 
a = 2.25 



233 



LENS FRONT VIEW 
OBJECT PLANE 



o 

3 

(O r- 
0^3 



Q_ 

(X 



(X 

r 



a: 

CO 



3 




J—J : ! LJ_ ! ^ I i A ^ L-i 



-I 



LO 


LxJ 


O 


d 




d 


O 


1 — 


O 

1 


CJ 


1 


LU 




d 




CD 


LO 


O 

I 


O 


1 

rvj 



o 



'_n 




j N ' I 1 - J ( 9 0 “ L 



2 34 



Figure E-26 . Grid Plane at a =0.3 for Lens of Figure E-25 
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Figure E-27. Spot Diagram for Grid of Figure E-26 
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E-28. Encircled Energy of Figure E-27 
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Figure E-30. Grid Plane at a =0.3 for Lens of Figure E-29 
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Figure E-32. Encircled Energy of Figure E-31 
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Figure E-33. GRIN Lens Shape at -10%, OB = 0.30, 
a = 2.25 
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Fiaure E-35. Spot Diagram for Grid of Figure E-34 
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Figure E-37. GRIN Lens Shape at -25%, OB = 0.35, 
a = 2.25 
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Figure E-38. Grid Plane at a =0.3 for Lens of Figure E-37 
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Figure E-39. Spot Diagram for Grid of Figure E-38 
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Figure E-40. Encircled Energy of Figure E-39 
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Figure E-42. Grid Plane at a =0.3 for Lens of Figure E-41 
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Figure E-43. Spot Diagram for Grid of Figure E-42 
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Figure E-44. Encircled Energy of Figure E-43 
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Figure E-45. GRIN Lens Shape at -10%, OB = 1.00 
a = 2.25 
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Figure E-4b. Grid Plane at a =0,3 for Lens of Figure 
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Figure E-47. Spot Diagram for Grid of Figure E-46 
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Figure E-48. Encircled Energy of Figure E-47 
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Figure E-51. Spot Diagram for Grid of Figure E-50 
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Figure E-52, Encircled Energy of Figure E-51 
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Figure E-53, GRIN Lens Shape at +5%, OB = 1.00, 
a = 2.25 
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Figure E-54. Grid Plane at a =0.3 for Lens of Figure E-53 
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Figure E-55. Spot Diagram for Grid of Figure E-54 
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Figure E-58. Grid Plane at a =0.3 for Lens of Figure 
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Figure E-59. Spot Diagram for Grid of Figure E-58 
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Figure E-62. Grid Plane at a =0.3 for Lens of Figure 
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Figure E-63. Spot Diagram for Grid of Figure E-62 
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Figure S-64. Encircled Energy of Figure E-63 
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Figure E-65. GRIN Lens Shape at +50%, OB - 1.00, 
a = 2.25 
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Figure E-66. Grid Plane at a = 0.3 for Lens of Figure 
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Figure E-67. Spot Diagram for Grid of Figure E-66 
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Figure E-70. Grid Plane at a =0.3 for Lens of Figure E-69 
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Figure E-72. Encircled Energy of Figure E-71 
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Figure E-75. Spot Diagram for Grid of Figure E-74 
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Figure E-76. Encircled Energy of Figure E-75 
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Figure E-77. GRIN Lens Shape at +5%, OB - 4.00, 
a = 2.25 
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Figure E-79. Spot Diagram for Grid of Figure E-78 
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Figure E-82. Grid Plane at a =0.3 for Lens of Figure 
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Figure E-83. Spot Diagram for Grid of Figure E-82 
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Figure E-85, GRIN Lens Shape at +25%, 
a = 2.25 
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Figure S-87. Spot Diagram for Grid of Figure E-86 



295 



SPOT DIflGRRM ENERGY DENSITY 












1 . 0 



0.8 



0. S 



0.4 



0.2 



0.0 




0.0 0.1 ].2 0.3 0.4 0.5 



RRD I uS 



u 






CENTROID 



Figure E-88, Encircled Energy of Figure E-87 
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Figure E-89 . GRIN Lens Shape at +50%, OB = 4.00, 
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Figure F-3. Spot Diagram for Grid of Figure F-2 
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Figure F-4. Encircled Energy of Figure F-3 
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Figure F-5. GRIN Lens Shape at -10%, OB = 0.05, 
a = 9.00 
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Figure F-8. Encircled Energy of Figure F-7 
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Figure F-9. GRIN Lens Shape at -25%, OB = 0.05, 
a = 9.00 
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Figure F-10 . Grid Plane at a =0.3 for Lens of Figure 
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Figure F-11. Spot Diagram for Grid of Figure F-10 



311 



I NfRGY OrNbl 1 1 



SPOT OlflGRflM ENERGY DENSITY 




npc:us -ROM :mpoe cem'^roid 



Figure F-12. Encircled Energy of Figure F-11 
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Figure F-13. GRIN Lens Shape at -50%, OB = 0.05, 
a = 9.00 
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Figure F-14. Grid Plane at a =0.3 for Lens of Figure F-13 
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Figure F-15. Spot Diagram for Grid of Figure F- 
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Figure F-17. GRIN Lens Shape at +5%, OB = 0.05, 
a = 9.00 
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Figure F-19. Spot Diagram for Grid of Figure F-18 
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Figure F-20. 
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Figure F-21, GRIN Lens Shape for +10%, OB = 0.05, 
a = 9.00 
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Figure F-22. Grid Plane at a =0.3 Radians for Lens of Figure F-21 
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Figure F-23. Spot Diagram for Grid of Figure F-22 
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Figure F-24. 



Encircled Energy of Figure F-23 
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Figure F-26 . Grid Plane at a = 0.3 for Lens of Figure 
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Figure F-27. Spot Diagram for Grid of Figure F-26 
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Figure F-29. GRIN Lens Shape at +50%, OB = 0.05, 
a = 9.00 
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Figure F-32. Encircled Energy of Figure F-31 
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Figure F-33. GRIN Lens Shape at -5%, OB = 0.20, 
a = 9.00 
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Figure F-35. Spot Diagram for Grid of Figure F-34 
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Figure F-36. Encircled Energy of Figure F-35 
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Figure F-37. GRIN Lens Shape at -10%, OB = 0.20, 
a = 9.00 
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Figure F-38. Grid Plane at a = 0 . 3 for Lens of Figure F-37 
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Figure F-40 . Encircled Energy of Figure F-39 
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Figure F-41. GRIN Lens Shape at -25%, OB = 0.20, 
a = 9.00 
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Spot Diagram for Grid of Figure F-42 
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Figure F-44. 



Encircled Energy of Figure F-43 
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Figure F-46. Grid Plane at a =0.3 for Lens of Figure F-45 
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Figure F-47. Spot Diagram for Grid of Figure F-46 
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Figure F-49. GRIN Lens Shape for + 5 %, OB = 0.20, 
a = 9.00 
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Figure F-50, Grid Plane at a = 0 . 3 for Lens of Figure F-49 
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Figure F-51. Spot Diagram for Grid of Figure F-50 
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Figure F-52. Encircled Energy of Figure F-51 
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Figure F-53. GRIN Lens Shape at +10%, OB = 0.20, 
a = 9.00 
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Figure F-54. Grid Plane at a =0.3 for Lens of Figure F-53 
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Figurs F~55 . Spot Diagirain for Grid of Figuro F~54 
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Figure F-56. Encircled Energy of Figure F-55 
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Figure F-5 7. GRIN Lens Shape at +25% r OB = 0.20, 
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Figure F-58. Grid Plane at a =0.3 for Lens of Figure 
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Figure F-59. Spot Diagram for Grid of Figure F-58 
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Figure F-60. Encircled Energy of Figure F-59 
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Figure F-62. Grid Plane at a =0.3 for Lens of Figure 
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Figure F-66. Grid Plane at a =0.3 for Lens of Figure F-65 
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Figure F-67. Spot Diagram for Grid of Figure F-66 
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Figure F-68. Encircled Energy of Figure F-67 
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Figure F-70. Grid Plane at a =0.3 for Lens of Figure F-69 
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Figure F-71. Spot Diagram for Grid of Figure F-70 
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Figure F-74. Grid Plane at a =0.3 for Lens of Figure 
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Figure F-75. Spot Diagram for Grid of Figure F-74 
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Figure F-79 . Spot Diagram for Grid of Figure F-78 
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Figure F-80. Encircled Energy of Figure F-79 
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Figure F-83. Spot Diagram for Grid of Figure F-82 
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Figure F-84. Encircled Energy of Figure F-83 
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Figure F-86. Grid Plane at a =0.3 for Lens of Figure 
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Figure F-87. Spot Diagram for Grid of Figure F-86 
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Figure F-88. Encircled Energy of Figure F-87 
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Figure F-91. Spot Diagram for Grid of Figure F-90 
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Figure F-93. GRIN Lens Shape at +50%, OB = 1.00 
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Figure F-94. Grid Plane at a =0.3 for Lens of Figure F-93 
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Figure F-95. Spot Diagram for Grid of Figure F-94 
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Figure F-97. GRIN Lens Shape at -5%, OB - 2.00, 
a = 9.00 
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Figure F-98. Grid Plane at a =0.3 for Lens of Figure F-97 
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Figure F-99 . Spot Diagram for Grid of Figure F-98 
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Figure F-100. Encircled Energy of Figure F-99 
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Figure F-102. Grid Plane at a =0.3 for Lens of Figure F-101 
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Figure F-103. Spot Diagram for Grid of Figure F-102 
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Figure F-105. GRIN Lens Shape at -25%, OB = 2.00, 

a = 9.00 
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Figure F-106. Grid Plane at a =0.3 for Lens of Figure F-105 
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Figure F-107. Spot Diagram for Grid of Figure F-106 
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F-108. Encircled Energy of Figure F-107 
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Figure F-109. GRIN Lens Shape at -50%, OB = 2.00, 

a = 9.00 
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Figure F-110. Grid Plane at a =0.3 for Grid of Figure F-109 
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Figure F-113. GRIN Lens Shape at +5%, OB = 4.00, 

a - 9.00 
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Figure F-114 . Grid Plane at a =0.3 for Lens of Figure F-113 
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-115. Spot Diagram for Grid of Figure F-114 
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Figure F-116. Encircled Energy of Figure F-115 
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Figure F-117. GRIN Lens Shape at +10%, OB - 4.00, 

a = 9.00 
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Figure F-118. Grid Plane at a =0.3 for Lens of Figure 
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Figure F-119. Spot Diagram for Grid of Figure F-118 
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Figure F-121. GRIN Lens Shape at +25%, OB = 4.00, 

a = 9.00 
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Figure F-122. Grid Plane at a =0.3 for Lens of Figure F-121 
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Figure F—123. 



Spot Diagram for Grid of Figure F-122 
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Figure F-124. Encircled Energy of Figure F-123 
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Figure F-125. GRIN Lens Shape at +50%, OB = 4.00, 

a = 9.00 
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Figure F-126. Grid Plane at a =0.3 for Lens of Figure F-125 
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Figure F-127. Spot Diagram for Grid of Figure F-126 
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Figure F-128. Encircled Energy of Figure F-127 
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Figure G-1. "Best" GRIN Lens Shape with 50% Gradient, 
OB = 0.05 and a = 9.00 in the F/2 
Configuration 
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Figure G-2. Grid Plane at a =0.0 for Lens of Figure 
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Figure G-4 . Grid Plane at a =0.1 for Lens of Figure 
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Figure G-5. Spot Diagram for Grid of Figure G-4 
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Figure G-6. Encircled Energy of Figure G-5 
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E'igure G-7. Grid Plane at a =0.2 for Lens of Figure 
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Figure G-8. Spot Diagram for Grid of Figure G-7 
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Figure G-9 . Encircled Energy of Figure G-8 
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Figure G-10 . Grid Plane at a =0.4 for Lens of Figure 



IMRGE PLRNr 
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Figure G-12. Encircled Energy of Figure G-11 
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Figure G-13. Grid Plane at a =0.5 for Lens of Figure 
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Figure G-14 . Spot Diagram for Grid of Figure G-13 
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Figure G-15. Encircled Energy of Figure G-14 
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Figure G-16. Grid Plane at a =0.7 for Lens of Figure 
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Figure G-17. Spot Diagram for Grid of Figure G-16 
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Figure G-18. Encircled Energy of Figure G-17 
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Figure G-19 . Grid Plane at a =0.7 for Lens of Figure 
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Figure G-20. Spot Diagram for Grid of Figure G-19 
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Figure G-21. Encircled Energy of Figure G-20 
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Figure G-22. Grid Plane at a = 0.8 for Lens of Figure 
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Figure G-23. Spot Diagram for Grid of Figure G-22 
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Figure G-24 . Encircled Energy of Figure G-23 
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Figure H-1. "Best" GRIN Lens Shape with 50% Gradient, 
OB = 0.05, and a = 9.00 in the F/1 
Configuration 
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Figure H-2. Grid Plane at a =0.0 for Lens of Figure 
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Figure H-3. Spot Diagram for Grid of Figure H-2 
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Figure H-4 . Grid Plane at a =0.1 for Lens of Figure 
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Figure H-6. Encircled Energy of Figure H-5 
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Figure H-7. Grid Plane at a =0,2 for Lens of Figure H- 
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Figure H-8. Spot Diagram for Grid of Figure H-7 
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Figure H-9 . Encircled Energy of Figure H-8 
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Figure H-10. Grid Plane at a =0.3 for Lens of Figure 
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Figure H-11. Spot Diagram for Grid of Figure H-10 
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Figure H-12. 



Encircled Energy of Figure H-11 



464 



LENS FRONT VIEN 
OBJECT PLANE 



O 

O 

O 

O 

O 

IT rvj 
O nn 
II 

It 

Q_ 

(X Ul 




O 

II 

cn 

>- 

cn 

GC 

Q 



j. ! ! I I I i i ! I j- ! .t ..L. > ! ! ! 1. 





O 

O 

I 




LO 



LD 

O 



O 

o 

I 



LD 

O 

I 



dMbnd iDdPdO Jv 



o 



I 



LiJ 

CE 

_1 

Q_ 



I— 

LJ 

LU 

~1 

CD 

O 

I 

rvj 



I 

a: 



465 



Figure H 13. Grid Plane at a = 0.4 for Lens of Figure 
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Figure H-14 . Spot Diagram for Grid of Figure H-13 
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Figure H-15. Encircled Energy of Figure H-14 
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Figure H-16. Grid Plane at a =0.5 for Lens of Figure 
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Figure H-17. Spot Diagram for Grid of Figure H-16 
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Figure H-18. Encircled Energy of Figure H-17 
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Figure H-19. Grid Plane at a =0.6 for Lens of Figure H- 
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Figure H-20. Spot Diagram for Grid of Figure H-19 
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Figure H-22. Grid Plane at a =0.7 for Lens of Figure 
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Figure H-23. Spot Diagram for Grid of Figure H-22 
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Figure H-24 . Encircled Energy of Figure H-23 
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Figure H-25. Grid Plane at a =0.8 for Lens of Figure 
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Figure H-26. Spot Diagram for Grid of Figure H-25 
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Figure H-27. Encircled Energy of Figure H-26 
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